Disordered weak and strong topological insulators 
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Robustness against disorder is a defining property of the topological quantum phenomena. Here, 
we highlight unexpected robustness of transport characteristics found in a lattice model of disordered 
three-dimensional Z2 topological insulator. We have studied numerically the global phase diagram 
of this model yielding both the weak and strong topological insulator (WTI and STI) phases to 
quantify how they collapse as a function of disorder. We have found that the average two-terminal 
conductance is quantized both in the bulk and slab geometries. This indicates that not only the 
surface Dirac cones in the STI and WTI phases but also bulk Dirac cones emergent at the phase 
boundaries exhibit robustness against disorder. We have also studied the Lyapunov exponents in 
the quasi one-dimensional geometry to verify that both the STI and WTI phases are indeed stable 
up to a finite strength of disorder. 

PACS numbers: 73.20.-r, 73.20.Fz, 71.23.-k, 71.30.+h 
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The Z2 topological insulator [1-3] is a time-reversal 
invariant counterpart of the quantum Hall state. The 
former realizes due to spin-orbit coupling preserving the 
time reversal symmetry (TRS), while the latter occurs 
under a strong magnetic field breaking the TRS. As the 
quantum Hall state is characterized by an integral (Z- 
type) topological invariant (Chern number), the Z2 topo- 
logical insulator is distinguished from "ordinary" insula- 
tors (01) by an analogous Z2 topological number, taking 
only two possible values: and 1. For three-dimensional 
(3D) Z2 topological insulators, there exist four of such 
Z2 topological invariants [4-6] implying the existence of 
the so-called strong and weak topological insulators (STI 
and WTI) [7]. Here, the term "strong" or "weak" is 
prescribed by the robustness of the state against per- 
turbations that preserve the system's underlying TRS. 
Though WTI is a priori considered weak, it is definitely 
an open issue both theoretically and experimentally to 
quantify how weak the WTI actually is against disorder 
[8, 9]. This should be naturally discussed in the light 
of how robust the STI is in the presence of disorder. In 
the following we address these issues by studying numeri- 
cally the global phase diagram of such a 3D disordered Z2 
topological insulator implemented on a lattice employing 
the transfer matrix method. 

Topological insulators are characterized by a bulk 
topological invariant encoding their nontrivial band 
structure, but through the bulk/ surface correspondence 
this is equivalent to say that they exhibit topologically 
protected gapless surface states. While an STI exhibits a 
single helical Dirac cone that is protected, a WTI man- 
ifests generally an even number of Dirac cones that can 
be null depending on the orientation of the surface [7] . In 
the study of graphene, the unusual robustness of Dirac 
electrons (especially in the case of a single Dirac cone) 
against disorder has been widely recognized. As a con- 
sequence of the absence of backward scattering [10], the 
Dirac electrons do not localize [11, 12]. However, in the 
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FIG. 1: (Color online) The "global phase diagram" of dis- 
ordered Z2 topological insulator in the (mo, W^)-plane deter- 
mined by the behavior of two-terminal conductance. Solid 
lines on the phase boundaries are guides to the eyes. Dotted 
lines indicates the value of parameters relevant in Figs. 2 and 
3. Metallic (M) phase lies in the intermediate range of disor- 
der strength, typically 10 < W < 25 in the parameter range 
of mo shown in this figure. 



presence of valleys (even number of Dirac cones) as in 
the case of graphene, they do localize mediated by inter- 
valley scatterings [13]. Does this mean that an STI con- 
tinues to be an STI for an arbitrary strength of short- 
ranged potential disorder that induces an inter-valley 
scattering, while a WTI simply collapses on switching on 
of the disorder? The reality is quite different from this 
naive speculation as shown in the phase diagram depicted 
in Fig. 1. This phase diagram established by a combina- 
tion of the study of the averaged two-terminal conduc- 
tance and of the quasi-lD decay length in the transfer 
matrix approach is intended to sketch the main results of 
the paper. In the phase diagram toq is the magnitude of 
the bulk energy gap at the F point (when mo is negative, 
the gap is inverted), while W is the strength of disorder. 
Notice also that in the model employed [see Eq. (1)], a 
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FIG. 2: Two-terminal conductance (G) (a) in the slab and 
(b) in the bulk geometries. (G) is plotted against mo for a 
system of size L = 10 under different strength of disorder: 
W = 1 (green), W = Z (red), and W = 5 (blue). In the bulk 
geometry PBC's are imposed in both the x- and {/-directions, 
while in the slab geometry, FBC is applied to the a::-direction. 



WTI phase with weak indices [7] (vQ^viViv-i) — (0,111) 
appears in the regime: —4 < mo < —2. The entire 
phase diagram is symmetric with respect to mo = —3. 
To identify the nature of different phases and the loca- 
tion of the phase boundaries in the (mo, Vl^)-plane, use 
of different geometries (i.e., bulk vs. slab) is shown to be 
crucial. While a plateau of the conductance in the slab 
geometry characterizes the nature of the corresponding 
phase [e.g., Fig. 2(a)], the phase boundaries are marked 
by a peak of the conductance in the bulk geometry [e.g.. 
Fig. 2(b)]. Under the breaking of translational invari- 
ance by disorder, standard techniques [7] for calculating 
the topological invariants fail. Yet, the above behaviors 
of the conductance clearly distinguish different topolog- 
ical phases, providing us with sufficient information for 
establishing the phase diagram depicted in Fig. 1. 

Let us begin by pointing out that both STI and WTI 
phases in the clean limit {W — 0) survive in the presence 
of finite disorder, but collapse into a metallic phase (the 
M-phase in Fig. 1) at a certain strength of disorder. We 
have confirmed the robustness of these insulating phases 
by studying the system size dependence of the average 
conductance (G) . The existence of surface helical Dirac 
cones and their nature are revealed by studying the gen- 
eralized quasi- ID decay length ^i. 

As for the global structure of the phase diagram, note 
that the STI exhibits a direct boundary with the 01 phase 
[14], without being intervened by an extension of the 



symplectic metal phase [15]. This is quite contrary to 
the case of the phase diagram for the 2D version of our 
model [16, 17]. In the weakly disordered regime, a couple 
of marked features are to be mentioned. First the WTI 
phase has an internal structure; it is divided into WTI 
and defeated WTI (DWTI) regions, reflecting the change 
of the system size dependence of the conductance. This 
indicates that a WTI phase is in a sense indeed weak 
compared with the STI. We leave detailed description of 
the DWTI region to the second half of the paper. An- 
other remarkable feature in this regime is the shape of 
the phase boundaries between different insulating phases. 
The positions of these phase boundaries are determined 
by the behavior of average conductance in the bulk and 
in the slab geometries (see Fig. 2). In the case of STI/OI 
boundary, a gradual shift of this position (to the 01 side) 
as increasing disorder represents a feature often referred 
to as a "topological Anderson insulator" [18-20]. A sim- 
ilar gradual shift of the phase boundary can be also seen 
at the WTI/STI boundary. On this side of the phase, 
the STI turns out to be less expansive, invaded by the 
WTI. This tendency of the phase boundary is consistent 
with the SCBA calculation [21]. 

Another unexpected feature we found in this study is 
the quantization of the peak conductance at the phase 
boundaries [see, e.g.. Fig. 2(a)]. At the phase bound- 
ary, the bulk energy gap closes and in the model stud- 
ied 3D Dirac cones emerge in the spectrum, dominating 
the transport characteristics of the system. The num- 
ber of such Dirac cones is 1 at the STI/OI boundary, 
while it is 3 at the WTI/STI boundary. As can be 
seen in Fig. 2(a), the height of the conductance peak 
is approximately given by twice this number in units of 
eV/i, i.e., G = 2 at the STI/OI boundary, and G = 6 
at the WTI/STI boundary. This is because a 3D bulk 
Dirac cone gives rise to two ID channels in the vicinity 
of i5 = 0. By studying the size dependence of G, we 
have verified that the height of the peak value of G ap- 
proaches indeed the quantized value, given by twice the 
number of Dirac cones. The quantization is an indication 
of a certain degree of robustness of the 3D Dirac cones 
against disorder. We mention that the fluctuation of G 
is strongly suppressed as approaching the peak. This is 
in contrast to the quantum (spin) Hall transition, where 
a finite universal fluctuation is observed [22] . 

As a concrete implementation of a 3D Z2 topological 
insulator on a lattice, we consider the following Wilson- 
Dirac type tight-binding Hamiltonian [23, 24], 
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(mo + 3r) ^4/3cx + ^t'x4l4Cx, (1) 



where c]^ and Cx are creation and annihilation operators 
on a site x, is a unit vector, ak and /3 are gamma 
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matrices in the Dirac representation 
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where au are PauU matrices and I2 is 2 x 2 identity ma- 
trix, t, r, mo are parameters, and disorder potential 
are uniformly distributed between —W/2 and W/2. We 
set r = 1, t = 2, as in Ref. [24]. This Hamiltonian belongs 
to the symmetry class All [25] (or, Dili for W = 0). 

The transfer matrix [26, 27] is given in terms of the 
wave function ij^n on a slice at z = n as 



H+ 



where Hn = {n\H\n) - E, H- = {n\H\n + 1) and H+ = 
(n + l\H\n). We set i? = in this study, though similar 
results are obtained for E = 0.05. 

To determine the phase boundaries between differ- 
ent insulating phases, we calculate the (average) two- 
terminal conductance, using the Landauer formula [22]. 
The transport between the left and right terminals is de- 
scribed in terms of the scattering matrix S defined as 



^out 



= S 



s = 



r t' 
t r' 



(3) 



where V'l(^^"'^ denotes the incoming (outgoing) state on 
the left (right) terminal, and t and t' (r and r') are 
transmission (reflection) matrices. The conductance G 
in units of e^/h is given by G = Tr{tH). To find the 
explicit form of the S^-matrix, we employ a transfer ma- 
trix for a system of length L. The actual computation 
has been done in a cubic geometry oi L x L x L sites 
with electrodes attached to the z-direction. We assume 
that the electrodes consist of perfectly-conducting ID 
wires as in network models [22], so that details of the 
electrodes do not affect the conductance. 

Typical examples of the calculated conductance curves 
are shown in Fig. 2. In the upper panel (a) the con- 
ductance is calculated in the "slab" geometry, i.e., fixed 
boundary condition (FBC) in the x-direction and peri- 
odic boundary condition (PBC) in y-direction (note z- 
direction is the direction of transport). In the lower panel 
(b), PBC's are imposed in both the x- and y-directions 
("bulk" geometry). The conductance is indeed very sen- 
sitive to the change of these boundary conditions. In the 
slab geometry, on the contrary, the conductance shows a 
plateau behavior, also quantized at (G) = 2 in the STI, 
and at (G) ~ 4 in the WTI phases. In the bulk geom- 
etry, the conductance exhibits a sharp peak feature at 
the phase boundary between different insulating phases, 
while inside the insulating phases, irrespective of their 
topological non-triviality (either in the STI, WTI, or 01 
phases), the conductance tends to vanish in the thermo- 
dynamic limit. Note that if one considers the robustness 
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FIG. 3: Two-terminal conductance as a function of disorder 
with mo = —2.5 for different system sizes: L = 10 (dotted), 
L — 14: (dashed), and L = 20 (sohd). The upper panel (a) is 
for the slab geometry. The weak topological surface states are 
stable for a certain value of disorder (about W < 4, which is 
actually not small), and defeated by disorder around W = 6.0. 
For W >7, the system enters metallic phase. The lower panel 
(b) shows the conductance in bulk geometry. There are no 
peak that should appear at the topological phase transition 
point around IF ~ 4. The statistical error is less than 0.02. 



of the helical surface Dirac cones (or, at least of a single 
Dirac cone) emergent on the surface of a slab, quanti- 
zation of these conductance plateaus is an expected be- 
havior up to a different degree of the quantization in the 
STI and WTI phases. In contrast, apparent quantization 
of the height of conductance peaks in the bulk geometry 
[Fig. 2(a)] was unexpected, since this quantization stems 
from the bulk Dirac cones. The positions of the conduc- 
tance peaks show little dependence on L and have been 
used to determine the phase boundaries of Fig. 1. 

To reinforce the validity of the phase diagram, we have 
also studied the Lyapunov exponents in the quasi-one- 
dimensional geometry. A Lyapunov exponent 7i is de- 
fined as 



7, = iim 



(4) 



where Ai is a positive eigenvalue of the matrix T^T, and 
T = Tm ■ ■ ■ T2T1 is a product of the M transfer matrices 
of dimension SL^Ly [26]. Due to the current conserva- 
tion, Xi occurs in reciprocal pairs and due to the Kramers 
degeneracy, they are doubly degenerate. Keeping this in 
mind, we arrange the exponents in the decreasing order, 

72L,L„ > 72L,L„-1 > ■ • • > 72 > 71 > > -71 > ■ • ■ > 
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FIG. 4: (Color online) Decay lengths (solid lines) and ^2 
(dashed lines) in the STI (mo,W) = (-1,1) (blue o), the 
WTI (mo,W) = (-2.5,3.5) (red □), and DWTI {mo,W) = 
( — 2.5,6) (green x) regions in the slab geometry with system 
width Ly — 6. The error bars are less than 5%. 

—l2L^Ly- The smallest positive Lyapunov exponent 71 
is identified as the quasi-lD decay length ^ by the cor- 
respondence, ^ = 1/71. Here, we extend this to higher 
Lyapunov exponents [28, 29]: 

= -. (5) 

li 

In the metallic phase, all the generalized decay length ^^'s 
uniformly increase as functions of the size of the cross sec- 
tion Lx and Ly. In the 01 phase, the decay lengths are 
insensitive to L^ and Ly. In the topological insulating 
phases with slab geometry, higher decay lengths behave 
similarly to those of the 01 phase, while the behavior of 
the largest few ^^'s is clearly distinguishable from those of 
the 01 phase. They are significantly large, but in contrast 
to the metallic states, they show specific dependence on 
the system's thickness L^. These behaviors reflect the na- 
ture of the 2D surface states. Indeed, the number of these 
special ^i's corresponds to the number of Dirac cones on 
a surface (1 for STI and 2 for WTI). 

Let us carefully see the transport properties of the STI 
and WTI phases (in Fig. 1) in the presence of disorder 
by using the conductance and decay length. In the STI 
phase, the conductance approaches the quantized value 
G = 2 irrespective of the disorder strength, and form a 
conductance plateau [see Fig. 2(a)]. In this plateau re- 
gion, the largest decay length ^1 increases rapidly (expo- 
nentially) with Lx (see, e.g.. Fig. 4), while the others ^i>2 
are of the same order of ^i's for ordinary insulators. This 
implies that there is a single robust Dirac cone on each 
surface and in the limit La, — )> cxo, it becomes perfectly 
conducting (c/. systems of an odd number of channels, 
see Ref. [30, 31]). We have noticed that the dependence 
of S^i on the system's width Ly is notably complicated. 
As Ly increases, higher order ^^'s start to exhibit a non- 
monotonic dependence. The number of such anomalous 
fi's increases with the increase of Ly in a certain interval 
of Ly, and each time this number changes by two per 
Dirac cone. Presumably, this is caused by opening of 



higher subbands of the helical Dirac cones as a ID chan- 
nel separated from the bulk bands. Thus in the ther- 
modynamic limit Lx,Ly,Lz — 00, the conductance may 
become infinite, instead of taking a quantized value. 

In the WTI phase, we have found that there are two 
regions where the conductance shows different depen- 
dence on the cube size L. For weak disorder {W < 4 
for Too — —2.5), the conductance increases and asymp- 
totically approaches a finite value G ~ 4 as L increases 
[see Fig. 3(a)]. In this weakly disordered regime, the de- 
cay lengths and ^2 tend to increase as ^1 in the STI 
phase, before being saturated as Lx increases (see Fig. 4) 
at a value ca. 50 times larger than ^3. The saturation 
of ^ is caused by a small but finite inter- valley scattering 
amplitude, which was absent in the STI phase. Com- 
bining these two observations, one can convince oneself 
that the surface states of a WTI indeed remain stable 
in this regime of a finite disorder strength. On the con- 
trary, once the disorder exceeds a certain value {W ~ 4), 
the conductance decays with increasing L (Fig. 3) and ^1 
and ^2 do not show an exponential rise. Such a behavior 
is indeed indistinguishable from the case of 01. What is 
then the nature of these two regimes, WTI and DWTI in 
Fig. 1? Are they distinct phases, and is there a transi- 
tion between the two? Our result suggests the following. 
The absence of a conductance peak between the WTI 
and DWTI regions [see Fig. 3(b)] implies that they are 
topologically identical. Actually, ^1 and ^2, which cor- 
respond to the surface states, are still large and distin- 
guishable from the bulk states (although they are smaller 
than those in the STI and WTI phases). In the DWTI 
region, surface conducting states are simply "defeated" 
by disorder. In this sense we name this region the de- 
feated WTI or DWTI region. The fate of the WTI and 
DWTI regions in thermodynamic limit is highly nontriv- 
ial. We leave more substantial presentation on this point 
to a forthcoming publication. 

We have seen in this paper that the surface states of 
a WTI exhibit unexpected robustness, but are defeated 
by a sufficient strength of disorder. It is also shown that 
a single helical Dirac cone on the surface of an STI is 
not unrivaled. In the lattice model studied, it collapses 
under a finite strength of disorder, just as in the case of 
WTI surface states. Disorder can also drive an STI into 
a WTI. We have proposed that topologically non-trivial 
phases in the presence of disorder can be characterized 
by the number of the (generalized) decay length that is 
significantly larger than the bulk decay length. Typically, 
this number is 1 for STI and 2 for WTI/DWTI. 
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